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Abstract 

Quantum computers will work by evolving a high tensor power of 
a small (e.g. two) dimensional Hilbert space by local gates, which can 
be implemented by applying a local Hamiltonian H for a time t. In 
contrast to this quantum engineering, the most abstract reaches of 
theoretical physics has spawned "topological models" having a finite 
dimensional internal state space with no natural tensor product struc- 
ture and in which the evolution of the state is discrete, H = 0. These 
are called topological quantum filed theories (TQFTs). These exotic 
physical systems are proved to be efficiently simulated on a quantum 
computer. The conclusion is two-fold: 

1. TQFTs cannot be used to define a model of computation stronger 
than the usual quantum model 

2. TQFTs provide a radically different way of looking at quantum 
computation. The rich mathematical structure of TQFTs might 
suggest a new quantum algorithm. 
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1 Introduction 

A topological quantum field theory (TQFT) is a mathematical abstraction, 
which codifies topological themes in conformal field theory and Chern-Simons 
theory. The strictly 2— dimensional part of a TQFT is called a topological 
modular functor (TMF) . It (essentially) assigns a finite dimensional complex 
Hilbert space to each surface S and to any (self)-diffeomorphism h 

of a surface a linear (auto)morphism V{h) : — >■ V^(S'). We restrict 
attention to unitary topological modular functors (UTMF) and show that 
a quantum computer can efficiently simulate transformations of any UTMF 
as a transformation on its computational state space. We should emphasize 
that both sides of our discussion are at present theoretical: the quantum 
computer which performs our simulation is also a mathematical abstraction 
- the quantum circuit model (QCM) [D][Y]. Very serious proposals exist for 
realizing this model, perhaps in silicon, e.g. [Ka], but we will not treat this 
aspect. 

There is a marked analogy between the development of the QCM from 
1982 Feynman [Fey] to the present, and the development of recursive func- 
tion theory in 1930's and 1940's. At the close of the earher period, "Church's 
thesis" proclaimed the uniqueness of all models of (classical) calculation: re- 
cursive function theory, Turing machine, A-calculus, etc.... The present paper 
can be viewed as supporting a similar status for QCM as the inherently quan- 
tum mechanical model of calculation. The modern reconsideration of com- 
putation is founded on the distinction between polynomial time and slower 
algorithms. Of course, all functions computed in the QCM can be computed 
classically, but probably not in comparable time. Assigning to an integer 
its factors, while polynomial time in QCM [Sh] is nearly exponential time 

1 /3 

~ e" according to the most refined classical algorithms. The origin of this 
paper is in thought [Fr] that since ordinary quantum mechanics appears to 
confer a substantial speed up over classical calculations, that some principle 
borrowed from the early, string, universe might go still further. Each TQFT 
is an instance of this question since their discrete topological nature lends 
itself to translation into computer science. We answer here in the negative by 
showing that for a unitary TQFT, the transformations V{h) have a hidden 
poly-local structure. Mathematically, V{h) can be reahzed as the restriction 
to an invariant subspace of a transformation Y[ 9% on the state space of a 
quantum computer where each Qi is a gate and the length of the composi- 
tion is linear in the length of /i as a word in the standard generators, "Dehn 
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twists" of the mapping class group = diffeomorphisms (S)/identity compo- 
nent. Thus, we add evidence to the unicity of the QCM. Several variants and 
antecedents of QCM, including quantum Turing machines, have previously 
be shown equivalent (with and without environmental errors) [Y]. 

From a physical standpoint, the QCM derives from Schrodinger's equa- 
tion as described by Feynman [Fey] and Lloyd [LI]. Let us introduce the 
model. Given a decision problem, the first or classical phase of the QCM is 
a classical program, which designs a quantum circuit to "solve" instances of 
the decision problem of length n. A quantum circuit is a composition U„ of 
operators or gates Qi G U(2) or U(4) taken from some fixed list of rapidly 
computable matrices^, e.g. having algebraic entries. The following short list 
suffices to efficiently approximate any other choice of gates [Ki]: 
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The gates are applied on some tensor power space (C^^^^C") of "k qubits" 
and models a local transformation on a system of k spin | particles. The 
gate g acts as the identity on all but one or two tensor factors where it acts 
as a matrix as above. This is the middle or quantum phase of the algorithm. 
The final phase is to perform a local von Neumann measurement on a final 
state V'finai = U„(?/'initiai) (or a commuting family of the same) to extract a 
probabilistic answer to the decision problem. (The initial states' 'j/'mitiai niust 
also be locally constructed). In this phase, we could declare that observing 
a certain eigenvalue with probability > | means "yes." We are interested 
only in the case where the classical phase of circuit design and the length of 
the designed circuit are both smaller than some polynomial in n. Decision 
problems which can be solved in this way are said to be in the computational 
class BQP: bounded-error quantum polynomial. The use of C^, the "qubit", 
is merely a convenience, any decomposition into factors of bounded dimension 
gives an equivalent theory. We say U is a quantum circuit over if all tensor 
factors have dimension = p. 

Following Lloyd [LI], note that if a finite dimensional quantum system, 
say (C^)®^, evolves by a Hamiltonian H, it is physically reasonable to assert 

^the i-th digit of each entry should be computable in poly{i) time 
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L 

that H is poly-local, H = ^i?^, where the sum has < poly{k) terms and each 

1=1 

Hi = Hi ® id, where Hi acts nontrivially only on a bounded number (often 
just two) qubits and as the identity on the remaining tensor factors. Now 
setting Plank's constant h = 1, the time evolution is given by Schrodinger's 
equation: = ^^'^'■tH -^vhereas gates can rapidly approximate [Ki] any local 
transformation of the form e^'^**^^ . Only the nonabelian nature of the unitary 

L 

group prevents us from approximating directly from as U e^'^*^*. However, 

1=1 

by the Trotter formula: 



where the error fl is measured in the operator norm. Thus, there is a good 
approximation to Ut as a product of gates: 

Vt = (e^-^^i . . . e^-^^^^)" + L^.o(^^y 

Because of the rapid approximation result of [Ki], in what follows, we will 
not discuss quantum circuits restricted to any small generating set as in the 
example above, rather we will permit a2x2or4x4 unitary matrix with 
algebraic number entries to appear as a gate. 

In contrast to the systems considered by Lloyd, the Hamiltonian in a 
topological theory vanishes identically, H = 0, a different argument - the 
substance of this paper - is needed to construct a simulation. The reader 
may wonder how a theory with vanishing H can exhibit nontrivial unitary 
transformations. The answer lies in the Feynman path-integral approach to 
QFT. When the theory is constructed from a Lagrangian (functional on the 
classical fields of the theory), which only involves first derivatives in time, 
the Legendre transform is identically zero [At], but may nevertheless have 
nontrivial global features as in the Aharonov-Bohm effect. 

Before defining the mathematical notions, we would make two comments. 
First, the converse to the theorem is an open question: Can some UTMF effi- 
ciently simulate a universal quantum computer? Fault tolerantly? We would 
conjecture the answer is yes to both questions. Second, we would like to 
suggest that the theorem may be viewed as a positive result for computa- 
tion. Modular functors, because of their rich mathematical structure, may 
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serve as higher order language for constructing a new quantum algorithm. In 
[Fr], it is observed that the transformations of UTMF's can readily produce 
state vectors whose coordinates are computationally difficult evaluations of 
the Jones and Tutte polynomials. The same is now known for the state vec- 
tor of a quantum computer, but the question of whether any useful part of 
this information can be made to survive the measurement phase of quantum 
computation is open. 

We would like to thank Greg Kupperberg and Kevin Walker for many 
stimulating discussions on the material presented here. 

2 Simulating Modular Functors 

We adopt the axiomatization of [Wa] or [T] to which we refer for details. 
Also see, Atiyah [At] , Segal [Se] , and Witten [Wi] . 

A surface is a compact oriented 2— manifold with parameterized bound- 
aries. Each boundary component has a label from a finite set 
C — {1, a, 6, c, . . . } with involution^, 1 = 1. In examples, labels might 
be representations of a quantum group up to a given level or positive energy 
representations of a loop group, or some other algebraic construct. Tech- 
nically, to avoid projective ambiguities each surface S is provided with a 
Lagrangian subspace L C and each diffeomorphism / : S — >^ E' is 

provided with an integer "framing/ signature" so the dynamics of the theory 
is actually given by a central extension of the mapping class group. Since 
these extended structures are irrelevant to our development, we suppress 
them from the notation. We use the letter i below to indicate a label set 
for all boundary components, or in some cases, those boundary components 
without a specified letter as label. 

Definition 1. A unitary topological modular functor (UTMF) is a functor 
V from the category of (labeled surfaces with fixed boundary parameteriza- 
tions, label preserving diffeomorphisms which commute with boundary param- 
eterizations) to (finite dimensional complex Hilbert spaces, unitary transfor- 
mations) which satisfies: 

1. Disjoint union axiom: V{Yi H F2, h n 4) = V{Yi, ii) O ¥{¥2, 4). 

2. Gluing axiom: let Yg arise from Y by gluing together a pair of boundary 
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circles with dual labels, x glues to x, then 

V{Y„i) = ^V{Y,{i,x,x). 

xeC 

3. Duality axiom: reversing the orientation of Y and applying ^ to labels 

corresponds to replacing V by V* . Evaluation must obey certain nat- 
urality conditions with respect to gluing and the action of the various 
mapping class groups. 

4. Empty surface axiom: 1^(0) = C 

5. Disk axiom: Va = V{D, a) = 




6. Annulus axiom: 14^;, = y{A, (a, 6)) 



!, if a = b 
0, if a^b 



7. Algebraic axiom: The basic data, the mapping class group actions and 
the maps F and S explained in the proof (from which V may be re- 
constructed if the Moore and Seiberg conditions are satisfied, see [MS] 
or [Wa] 6.4, I-I4) is algebraic over Q for some bases in of Va, K,2; 
and Vabc, where Vahc denote ^(P, (a, 6, c)) for a (compact) 3-punctured 
sphere P. 3-punctured spheres are also called pants. 



Comments : 

1) Prom the gluing axiom, V may be extended via dissection from simple 
pieces D, A, and P to general surfaces S. But must be canonically 
defined: this looks quite difficult to arrange and it is remarkable that 
any nontrivial examples of UTMFs exist. 

2) The algebraic axiom is usually omitted, but holds for all known exam- 
ples. We include it to avoid trivialities such as a UTMF where action 
by, say, a boundary twist is multiplication by a real number whose bi- 
nary expansion encodes a difficult or even uncomputable function: e.g. 
the i^^ bit is iff the i^^ Turing machine halts. If there are nontrivial 
parameter families of UTMF's, such nonsensical examples must arise 
- although they could not be algebraically specified. In the context of 
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bounded accuracy for the operation of diffeomorphisms V{h), axiom 7 
may be dropped (and simulation by bounded accuracy quantum cir- 
cuits still obtained), but we prefer to work in the exact context since 
in a purely topological theory exactness is not implausible. 

3) Axiom 2 will be particularly important in the context of a pants 
decomposition of a surface E. This is a division of E into a collection of 
compact surfaces P having the topology of 3— punctured spheres and 
meeting only in their boundary components which we call "cuffs." 

Definition 2. A quantum circuit U : (C^)®*^ (^^p^^fc ^. g^j^g 
simulate on W (exactly) a unitary transformation r : S S if there is a C- 
linear imbedding i : S C. (C^)®'^ invariant under U so that U oi = ior. The 
imbedding is said to intertwine r and U. We also require that i be computable 
on a basis in poly(k) time. 

Since we prove efficient simulation of the topological dynamics for UTMFs 
y, it is redundant to dwell on "measurement" within V, but to complete the 
computational model, we can posit von Neumann type measurement with 
respect to any efficiently computable frame in Vabc- The space above, 
later denoted X = C^, is defined hy X := © Vabc and the computational 

(a,6,c)6£3 

space W := X®*^. We have set S := V^(E) and assumed E is divided into 
k "pants," i.e. Euler class (E) = —k. Any frame JF extends to a frame 
for V^(E) via the gluing axiom once a pants decomposition of E is specified. 
Thus, measurement in V becomes a restriction of measurement in W. It may 
be physically more natural to restrict the allowable measurements on ^ (E) to 
cutting along a simple closed curve 7 and measuring the label which appears. 
Mathematically, this amounts to transforming to a pants decomposition with 
7 as one of its decomposition or "cufF' curves and then positing a Hermitian 
operator with eigenspaces equal to the summands of T^(E) corresponding 
under the gluing axiom to labels x on 7" and x on 7''", xeC 

A labeled surface (E, i) determines a mapping class group Ml = A^(E, i) — 
"isotopy classes of orientation preserving diffeomorphisms of E preserving la- 
bels and commuting with boundary parameterization." For example, in the 
case of an n-punctured sphere with all labels equal (distinct), M — SFB(n), 
the spherical framed braid group (At = PSFB(n), the pure spherical framed 
braid group) . To prove the theorem below, we will need to describe a gen- 
erating set S for the various A4's and within S chains of elementary moves 
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which will allow us to prepare to apply any S2 & S subsequent to having 
applied si G S. 

Each Ai is generated by Dehn-twists and braid-moves (See [B]). A Dehn- 
twist is specified by drawing a simple closed curve (s.c.c.) 7 on S, cutting 
along 7, twisting 2n to the right along 7 and then regluing. A braid-move Bs 
will occur only when a s.c.c. 6 cobounds a pair of pants with two boundary 
components of S: If the labels of the boundary components are equal then 
Bs braids them by a right vr- twist. In the case that all labels equal, there is 
a rather short list of D and B generators indicated in Figure 1 below. Also 
sketched in Figure 1 is a pants decomposition of diameter = (O log 61(E)), 
meaning the graph dual to the pants decomposition has diameter order log 
the first Betti number of S. 




Figure 1. 

The s.c.c. 7( S) label Dehn (braid) generators (and B^). Figure 
1 contains a punctured annulus A; note that the composition of oppositely 
oriented Dehn twists along the two "long" components of dA, 7 and 7' yield 
a diffeomorphism which moves the punctures about the loop 7. The figure 
implicitly contains such an A for each (7,p), where p is a preferred puncture. 
The 7 curves come in three types: 

(1) The loops at the top of the handles which are curves ("cuffs") of the 
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pants decomposition, 

(2) loops dual to type 1, and 

(3) loops running under adjacent pairs of handles (which cut through up 
to (9(log(6iS)) many cuffs). (See Figure 1, where cuffs are marked by 
a "c".) 

Each punctured annulus A is determined as a neighborhood (of a s.c.c. 7 
union an arc rj from 7 to p). To achieve general motions of p around S, we 
require these arcs to be "standard" so that for each p, 7ri(S,p) is generated 
by ■ 7 ■ ^7"^}, where S = S with punctures filled by disks, and the disk 
corresponding to p serving as a base point. This list of generators is only 
linear in the first Betti number of S. 

In the presence of distinct labels, many of the Bs are illegal (they permute 
unequal labels). In this case, quadratically many generators are required. 
Figure 2 displays the replacements for the -B's, and additional A's and D's . 




Figure 2. 

Figure 2 shows a collection of B's sufficient to effect arbitrary braiding 
within each commonly-labeled subset of punctures, a quadratically large col- 
lections of new Dehn curves {e} allowing a full twist between any pair of 
distinctly labeled punctures, (If the punctures are arranged along a convex 
arc of the Euclidean cell in S, then each e will be the boundary of a narrow 
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neighborhood of the straight hne segment joining pairs of dissimilarly labeled 

punctures.) and finally a collection of punctured annuli, which enable one 
puncture Pi from each label - constant subset to be carried around each free 
homotopy class from {7} (respecting the previous generation condition for 

Thus for distinct labels the generating sets are built from curves of type 
7,7',e and d by Dehn twists around 7,7', and e, braid moves around 5. 
Denote by cu, any such curve: cu E = {{7} U {7'} U {e} U {5}}. 

Since various a;'s intersect, it is not possible to realize all uj simultaneously 
as cuffs in a pants decomposition. However, we can start with the "base 
point" pants decomposition V indicated in Figure 1 (note 7 of type(l) are 
cuffs in V, but 7 of types (2) or (3) are not) and for any u find a short path 
of elementary moves: F and S (defined below) to a pants decomposition T>^ 
containing as a cuff. 

Lemma 2.1. Assume S 7^ S*^, disk, or annulus, and T> the standard pants 
decomposition sketched in Figure 1. Any uo as above, can he deformed through 
0(log6i(S)) F and S moves to a pants decomposition V^^ in which uj is a 
cuff. 

We postpone the proof of the lemma and the definition of its terms until 
we are partly into the proof of the theorem and have some experience passing 
between pants decompositions. 

Theorem 2.2. Suppose V is a UTMF and h : T, ^ T, is a diffeomorphism of 
length n in the standard generators for the mapping class group ofE described 
above (See Figures 1 and 2). Then there are constants depending only on V , 
c — c{y) andp — p{V) such that V{h) : V{T,) V{T,) is simulated (exactly) 
by a quantum circuit operating on '^qupits" of length < c ■ n ■ log 61(E). 

The collection {cuffs} refers to the circles along which the pants decom- 
position decomposes; the "seams" are additional arcs, three per pant which 
cut the pant into two hexagons. Technically, we will need each pant in V to 
be parameterized by a fixed 3-punctured sphere so these seams are part of 
the data in V; For simplicity, we choose seams to minimize the number of 
intersections with {cu}. 

The theorem may be extended to cover a more general form of input. The 
original algorithm [L] which writes a D„, (y, db S. C.C.J clS db word in standard 
generators is super-exponential. We define the combinatorial length of a, 
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i{a), to be the minimum number of intersections as we vary a by isotopy of 
a with {cuffs} U {seams}. The best upper-bound (known to the authors) to 
the length L of T>a as a word in the mapping class group spanned by a fixed 
generating set is of the form < super-exponential function f{i). For 

this reason, we consider as input V{h), where h is a composition of k Dehn 
twists on ai, . . .Ok and j braid moves along (3i, . . . j3j in any order. Then 
V{h) is costed as the sum of the combinatorial length of the simple closed 
curves needed to write h as Dehn twists and braid moves within the mapping 

j k 

class group, i{h) := S + S i{a,i). We obtain the following extension of 

i=l i=l 

the theorem. 

Extension: The map h^, : V{I1) — > V{Tj) is exactly simulated by a 
quantum circuit QC with length (QC) < lli{h) composed of algebraic 1 and 
2— qupit gates. 

Pre-Proof: Some physical comments will motivate the proof. V^(S) are 
quantized gauge fields on S (with a boundary condition given by labels i) 
and can be regarded as a finite dimensional space of internal symmetries. 
This is most clear when genus (S) = , S is a punctured sphere, the labeled 
punctures are "anyons" [Wil] and the relevant mapping class group is the 
braid group which moves the punctures around the surface of the sphere. An 
internal state tpeV{T,) is transformed to U(6)^/' G V{Ti) under the functorial 
representation of the braid group. For U(6) to be defined the braiding must 
be "complete" in the sense that the punctures (anyons) must return setwise 
to their initial position. Infinitesimally, the braiding defines a Hamiltonian H 
on VCE) ®E where E is an infinite dimensional Hilbert space which encodes 
the position of the anyons. The projection of H into V(S) vanishes which is 
consistent with the general covariance of topological theories. Nevertheless, 
when the braid is complete, the evolution U of if will leave V{Tj) invariant 
and it is U|y(s) = U which we will simulate. Anyons inherently reflect non- 
local entanglement so it is not to be expected that V{T?) has any (natural) 
tensor decomposition and none are observed in interesting examples. Thus, 
simulation of U as an invariant subspace of a tensor product (C^)®'= is the 
best result we can expect. The mathematical proof will loosely follow the 

•^Lee Mosher has informed us that the existence of linear bound f{l) (but without 
control of the constants) follows at least for closed and single punctured surfaces from his 
two papers [Ml] and [M2]. 
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physical intuition of evolution in a super-space by defining, in the braid case 
(identical labels and genus = 0), two distinct imbeddings "odd" and "even," 

odd 

\/(S)<^^(C^)® = W and constructing the local evolution by gates acting on 
the target space. The imbeddings are named for the fact that in the usual 
presentation of the braid group, the odd (even) numbered generators can be 
implemented by restricting an action on W to image odd V (T,) (^evenV {'E)). 

Proof: The case genus (S) = with all boundary components carrying 
identical labels (this contains the classical, uncolored Jones polynomial case 
[J] [Wi]) is treated first. For any number q of punctures [q = 10 in the illus- 
tration) there are two systematic ways of dividing S into pants (3— punctured 

spheres) along curves a = {ai, . . . , aqs} or along /3 = . . . , Pqs} so that 

a sequence of q F moves (6j— moves in physics notation) transforms a to (3. 




Figure 3. 
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Let -'^ = ©(a (, c)e£3 ^abc be the orthogonal sum of all sectors of the pants 
Hilbert space. Distributing over the tensor power X'^^'^~'^^ := 1^ is 
the sum over all labelings of the Hilbert space for IJ(5' — 2) pants. Choosing 
parameterizations, W is identified with both the label sum space (S^^^^) and 
sum (S ^). Now S is assembled from the disjoint union by gluing along a 

cut (3 

or (3 so the gluing axiom defines imbeddings i{a) and of ¥{12, i) as a 
direct summand of X^^''^'^^ = W . 

Consider the action of braid move about a. This acts algebraically as 
9{ai) on a single X factor of W and as the identity on other factors. This 
action leaves i{a) (y(E,^)) invariant and can be thought of as a "qupit" 
gate: 

e{ai) = ^(braidaj : X ^ X 

where dimension dim{X) — p. Similarly the action of ^(braid/j.) leaves i{(5) 

invariant. It is well known [B] that the union of loops a\J (5 determines a 
complete set of generators of the braid group. The general element which 
we must simulate by an action on is a word in braid moves on a's and 
/9's. Part of the basic data - implied by the gluing axiom for a UTMF is a 
fixed identification between elementary gluings: 



Fabcd ■ ® Ka6 ® Vxcd ' ^ Vybc V^da 

xeC yeC 

corresponding to the following two decompositions of the 4— punctured sphere 
into two pairs of pants (The dotted lines are pant "seams" , the uncircled num- 
ber indicate boundary components, the letters label boundary components, 
and the circled numbers order the pairs of pants.): 
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Figure 4. 

For each F, we choose an extension to a unitary map F' : X ^ X — > 
X ^X. Then extend F' to F by tensoring with identity on the g — 4 factors 
unaffected by F. The composition of q F's, extended to q F's, corresponding 
to the q moves illustrated in the case g = 10 by Figure 3. (For g > 10 imag- 
ine the drawings in figure 3 extended periodically.) These define a unitary 

transformation T : W W with T o i[a) = i{j3). The word u in the braid 
group can be simulated by r on W, where r is written as a composition of 
the unitary maps T, T~^, 6{ai), and 0{f3j). For example, 

/?5ai/?2'^«l«3 

would be simulated as 

T = T-^o e{l3^) oTo e{ai) o o e{f32^) oTo e{ai) o ^(ag). 

As described r has length < 2q length u. The dependence on q can be 
removed by dividing S into | overlapping pieces Sj, each Sj a union of 6 

consecutive pants. Every loop of a U /3 is contained well within some piece 
Ej so instead of moving between two fixed subspaces iaiV) and i/siV) C W, 
when we encounter a do constantly many F operations to find a new 
pants decomposition modified locally to contain Pj. Then 9{Pj) may be ap- 
plied and the F operations reversed to return to the a pants decomposition. 
The resulting simulation can be made to satisfy length r < 7 length u. This 
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completes the braid case with all bounding labels equal - an important case 
corresponding to the classical Jones polynomial [J]. 

Proof of Lemma: We have described the F— move on the 4— punc- 
tured sphere both geometrically and under the functor. The 5*— move is 
between two pants decompositions on the punctured torus T~ (Filling in the 
puncture, a variant of S may act between two distinct annular decomposition 
of T^. We suppress this case since, without topological parameter, there can 
be no computational complexity discussion over a single surface). 




By [Li] or [HT] that one may move between any two pants decomposi- 
tions via a finite sequence of moves of three types: F, S, and diffeomorphism 
M supported on the interior of a single pair of pants (appendix [HT]). To 
pass from V, our "base point" decomposition, to P^^, F and S moves alone 
suffice and the logarithmic count is a consequence of the log depth nest of 
cuff loops of V on the planar surface obtained by cutting E along type (1) 7 
curves. Below we draw examples of short paths of F and 5* moves taking T) 
to a particular Vi^. 

The logarithmic count is based on the proposition. 



Proposition 2.3. Let K be a trivalent tree of diameter = d and f he a 
move, which locally replaces 

{M} 

and with , then any two leaves of K 
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can be made adjacent by < d moves of type f. 

Passing from K to a. punctured sphere obtained by imbedding {K, uni- 
valent vertices) into (^i?^, i?^), thickening and deleting the boundary R^, the 
/ move induces the previously defined F move.DD 

Some example of paths of F, S moves: 




Figure 6. 
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Continuation of the proof of the theorem: For the general case, we 
compute on numerous imbeddings of V^(S) into W (rather than on two: 
ia{V{T,)^ and as in the braid case). Each imbedding is determined 

by a pants decomposition and the imbedding changes (in principle) via the 
lemma every time we come to a new literal of the word cu. Recall that cu & M., 
the mapping class group, is now written as a word in the letters (and their 
inverses) of type D^, D'^, D^, and Bs- Pick as a home base a fixed pants 
decomposition Vq corresponding to io(y{T,)^ C W. If the first literal is a 
twist or braid along the s.c.c. lu, then apply the lemma to pass through a 
sequence of F and S moves from Dq to Pi containing a; as a "cufP' curve. 
As in the braid case, choose extensions F and S to unitary automorphisms 
of W and applying V to the composition gives a transformation Ti of W 
such that ii = Ti o ig, being the inclusion V{T,) — > W associated with Vi. 
Now execute the first literal uui of as a transformation 0{ui), which leaves 
ii (V^(S)) invariant and satisfies: 9{cUi) oii = i^o V{uJi). Finally apply Tf ^ to 
return to the base inclusion io(y{T,)). The previous three steps can now be 
repeated for the second literal of uj: follow T^^o9{uJi)oTi by T2^o9{uj2)°T2. 
Continuing in this way, we construct a composition r which simulates uj on 
W: 

From lemma 2.1 the length of this simulation by one (corresponding to S 
and 9{LJi)^ and two (corresponding to F moves) qupit gates is proportional 
to n =length uj and log 61 (S), where p = dim(X). 

Proof of Extension: What is at issue is the number of preparatory moves 
to change the base point decomposition V to T>^ containing 7 = aj or f3i as 
a cuff curve 1 < i < k or j. We have defined the F and S moves rigidly, 
i.e. with specified action on the seams. This was necessary to induce a well 
defined action on the functor V. Because of this rigid choice, we must add 
one more move - an M move - to have a complete set of moves capable of 
moving between any two pants decompositions of a surface (compare [HT]). 
The M move is simply a Dehn twist supported in a pair of pants of the 
current pants decomposition; it moves the seams (compare chapter 5 [Wa]). 
Note that if M is a +1 Dehn twisit in a s.c.c. uj then, under the functor, 
V{M) is a restriction of 9{u;) in the notation above. 

As in [HT], the cuff curves of V may be regarded as level curves of a Morse 
function / : E ^ constant on boundary components which we assume to 
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have minimum complexity (= total number of critial points) satisfying this 
constraint. Isotope a (we drop the index) on E to have the smallest number 
of local maximums with respect to / and is disjoint from critical points of / 
on E. 

Now generically deform / in a thin annular neighborhood of 7 so that 7 
becomes a level curve. Consider the graphic G of the deformation ft, <t < 
1. For regular t the Morse function ft determines a pants decomposition: let 
the 1- complex K consist of E/ ~ where ,t ~ y if x and y belong to the same 
component of a level set of ft, and let L G K be the smallest complex to 
which K collapses relative to endpoints associated to boundary components. 
For example in figure 8, the top tree does not collapse at all while in the 
lower two trees the edge whose end is labeled, "local max" is collapsed away. 
The preimage of one point from each intrinsic 1— cell of L not containing a 
boundary point constitutes a {cuffs} determining a pants decomposition Vt. 
For singular to, let T>tQ-e and "Dto+e may differ or may agree up to isotopy. 
The only change in V occurs when i is a crossing point for index= 1 handles 
where the two critical points are on the same connected component of a 
level set f^^ij-). There are essentially only three possible "Cerf-transitions" 
and they are expressible as a product of 1, 2, or 3 F and S moves together 
with braid moves whose number we will later bound from above. The Cerf 
transitions on T> are shown in Figure 7, together with their representation as 
compositions of elementary moves. 
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ooo (Xg:) 

double critical level double critical level 




moves: F o M o F 
(as shown below) 



double critical level 




seams to 
standard position 

Figure 7. 

Critical points of become critical points of ft of the same index once 
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the deformation as passed an initial eo > 0, and before any saddle-crossings 
have occurred. Let P be a pant from the composition induced by / and S C 
7 n P an arc. Applying the connectivity criterion of the previous paragraph, 
we can see that flattening a local maxima can effect at most the two cuff 
circles which 6 meets, and these by elementary Cerf transition shown in 
Figure 8. 




Figure 8. Pulling 7 down yields F o F 



If 7 crosses the seam arcs then the transitions are of the Cerf type, pre- 
composed with M— moves to remove these crossings as shown in Figure 9. 
Dynamically seam crossings by 7 produce saddle connections in the Cerf 
diagram. 




pull y down, 
cancel local max ^ 

yields FoMo F 




Figure 9. 
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The total number of these twists is bounded by length (7). The number 
of flattening moves as above is less than or equal I7 fl cuffs | < lcngth(7). 
The factor of 11 in the statement allows up to 5 F, S, and M moves for 
expressing each Cerf singularity which arises in passing from Do to and 
the same factor of 5 to pass back from T>j to T>o again, while saving at least 
one step to implement the twist or build move along 7. This completes the 
proof of the extension. □ 

We should emphasize that although, we have adopted an "exact" model for 
the operation of the UTMF, faithful simulation as derived above does not 
depend on a perfectly accurate quantum circuit. Several authors have proved 
a threshold theorem [Ki], [AB], and [KLZ]: If the rate of large errors acting 
on computational qubits (or qupits) is small enough, the size of ubiquitous 
error small enough, and both are uncorrelated, then such a computational 
space may be made to simulate with probability > | an exact quantum cir- 
cuit of length = L. The simulating circuit must exceed the exact circuit in 
both number of qubits and number of operations by a multiplicative factor 
< poly (logL). 

3 Simulating TQFT's 

We conclude with a discussion about the three dimensional extension, the 
TQFT of a UTMF. In all known examples of TMF's there is an extension 
to a TQFT meaning that it is possible to assign a linear map 
V{H') subject to several axioms [Wa] and [T] whenever S and S' cobounds 
a bordism b (with some additional structure). The case of bordisms with 
a product structures is essentially the TMF part of the theory. Unitarity 
is extended to mean that if the orientation of the bordism b is reversed to 
6, we have bl — (6)*. It is known that a TMF has at most one extension 
to a TQFT and conjectured that this extension always exists. Non-product 
bordisms correspond to some loss of information of the state. This can be 
understood by factoring the bordism into pieces consisting of a product union 
a 2-handle: ExlUh. The 2-handle h has the form {D^ x I, dD^ x I) and 
is attached along the subspace dD'^ x I. The effect of attaching the handle 
will be to "pinch" off an essential loop a; on E and so replace an annular 
neighborhood of uj by two disks turning S into a simpler surface S'. It is 
an elementary consequence of the axioms that if6 = Sx/U/i then 6* is 
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a projector as follows: Let X> be a pants decomposition containing a; as a 
dissection curve. There are two cases: 

(1) CO appears as the first and second boundary components of a single pant 
called Pq or 

(2) CO appears as the first boundary component on two distinct pants called 
Pi and Pa- 



y(E) = 

= f (0 Vaac) (g) l^(E\Po, with label c on dsPo) ] , case (1), 
or 

= ( Kbc Vade J ( S\ (Pi U P2) , appropriate labels ) , case (2) . 

labels \aeC / \ / 

In case (2), there may be relation b = c and/or d = e depending on 
the topology of V. The map 6* is obtained by extending linearly from the 
projections onto summands: 



canonically 

0Kac Vul ^ Vi (easel) 

a,cC 

or 

. canonically _^ 

K6c(g)T4.e ^i,s(8)^idd = V,h(S^Vdd (case 2) 

If the orientation on b is reversed the unitarity condition implies that b 
determines an injection onto a summand with a formula dual to the above. 
Thus, any bordism's morphism can be systematically calculated. 

In quantum computation, as shown in [Ki], a projector corresponds to 
an intermediate binary measurement within the quantum phase of the com- 
putation, one outcome of which leads to cessation the other continuation of 
the quantum circuits operation. Call such a probabilistically abortive com- 
putation a partial computation on a partial quantum circuit. Formally, if we 
write the identity as a sum of two projectors: idy = IIq + Hi, and let Uq 
and Ui be unitary operators on an ancillary space A with Uo(|0)) = |0) and 
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Ui |0) = 1 1) . The unitary operator Ho ® Uq + Hi (g) Ui onV^A when apphed 
to \v) ® |0) is lllof) ® |0) + \Iliv) C?) |1) so continuing the computation only 
if the indicator |0) G A is observed simulates the projection Ho. 

It is clear that the proof of the theorem can be modified to simulate 2- 
handle attachments as well as Dehn twists and braid moves along s.c.c.'s lj 
to yield: 

Scholium 3.1. Suppose b is an oriented bordism from Sq to Hi, where Sj 
is endowed with a pants decomposition T>i. Let complexity (b) be the total 
number of moves of four types: F , S, M , and attachment of a 2-handle to a 
dissection curve of a current pants decomposition that are necessary to recon- 
struct b from (EqjVq) to (Si, Pi). Then there is a constant d{V) depending 
on the choice of UTQFT and p{V) as before (for the TQFTs underlining 
TMF) so that b^ : V{J1q) V{T,i) is simulated (exactly) by a partial quan- 
tum circuit over £P of length < d complexity (b) . 
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